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Abstract 

We calculate conductance through a quantum dot weakly coupled to metallic contacts by means of Keldysh out of 
equilibrium formalism. We model the quantum dot with the SU(2) Anderson model and consider the limit of infinite 
Coulomb repulsion. We solve the interacting system with the numerical diagrammatic Non-Crossing Approximation 
(NCA). We calculate the conductance as a function of temperature and gate voltage, from differential conductance 
(dl/dV) curves. We discuss these results in comparison with those from the linear response approach which can 
be performed directly in equilibrium conditions. Comparison shows that out of equilibrium results are in good 
agreement with the ones from linear response supporting reliability to the method employed. The discussion becomes 
relevant when dealing with general transport models through interacting regions. We also analyze the evolution of 
the curve of conductance vs gate voltage with temperature. While at high temperatures the conductance is peaked 
when the Fermi energy coincides with the energy of the localized level, it presents a plateau for low temperatures 
as a consequence of Kondo effect. We discuss different ways to determine Kondo's temperature and compare the 
values obtained in and out of equilibrium. 
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1. Introduction 

Since the first observation of the Kondo effect in 
semiconducting quantum dots (QD)[1] the study of 
transport through nanoscopic devices has inspired a 
rich variety of experimental and theoretical works. 
Nowadays, measurements of transport properties in 
such systems, such as current versus bias voltage and 
conductance, are the main focus of the experiments 
due to the interesting and unusual features observed 
[2,3]. 

The behavior of the conductance at different tem- 
peratures and for different gate voltages has been stud- 
ied in very general systems, including those showing 
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strong correlations. While at equilibrium almost exact 
numerical methods have been developed for the the- 
oretical treatment of this problems (numerical renor- 
malization group (NRG) [4] or exact diagonalization 
(ED)[5]), the ones for non-equilibrium conditions are 
still in progress. Among them, the Scattering Bethe 
Ansatz (SBA) [6] and the Time dependent Density Ma- 
trix Renormalization Group (t-DMRG) [7] are promis- 
ing. 

In this work we study the transport properties of an 
interacting QD using the Non-Crossing Approximation 
(NCA) in its non-equilibrium[8] and equilibrium [9] ver- 
sions. We consider mandatory the comparison beetwen 
both schemes to support reliability to the more general 
procedure dealing with an out of equilibrium calcula- 
tion. We discuss the results for conductance as function 
of bias and gate voltage, and moreover, the dependence 
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of transport properties with temperature. 



2. Model 

We study the transport properties through a quan- 
tum dot weak-coupled to metallic contacts and de- 
scribe the system with the Anderson model, 

H = ^2 ^k<T V c\ au c kav + ^2 Eddtdv + U d\dfd[di 

her is <r 

+ J2 VkvAvJ* + H - c -i (!) 

kva 

where Cfc CTI y(cJ, CTj/ ) is the destruction (creation) operator 
of an electron with momentum k, spin a and lead v = 
L (left) or R (right), and d a (dl) destroys (creates) an 
electron in the quantum dot. 

The non-interacting conduction electrons in the 
leads are treated as being in thermal and chemical 
equilibrium with their reservoirs, thus etav = tka — fJ-v, 
allowing for different chemical potentials p v in each of 
them. For the central region, we consider a spin de- 
generate localized level with energy Ed and Coulomb 
repulsion U. The leads and the dot are connected by 
means of hybridizations Vkau- 

The physical quantity accessible in transport mea- 
surements is the current. As shown by Meir and 
Wingreen[10], the current through a system described 
by the Hamiltonian Eq.(l) is given by 

+ nr Tr I MrL ~ rR)iG * {uj)/27V > (2) 

where F^ = 27r^ fc \Vkav\ 2 S(u — ek a v) is the hybridiza- 
tion function and is the Fermi function for the 
conduction electrons of the left(right) lead. The func- 
tions pd(u}) and (w) represent the spectral density 
and the lesser Green function of the central region 
respectively. The calculation of such Green functions 
must be done in presence of the leads and is a non- 
equilibrium problem which might be handled within 
Keldysh formalism. 

From the results of different applied bias voltages, 
differential conductance dl /dV curves can be obtained 
by differentiation of Eq. (2). The value at zero bias 
dI/dV\o — G(T) represent (the usual equilibrium con- 
ductance) at a given temperature. A simplified ana- 
lytical expression for G(T) can be obtained under the 
condition of proportional couplings, — aF„ [8] 

G (T) = 4^ J dcu (-?£\ T(uj)pd(uj), (3) 
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Fig. 1. QD spectral density for different applied bias voltages. 

where r ff = r£r^/(r£ +r^) is the effective hybridiza- 
tion and Pd(u) is calculated in equilibrium. Thus, in 
contrast to Eq. (2) ony equilibrium quantities enter Eq. 
(3). 



3. Results 

For the treatment of the model in and out of 
equilibrium, we use the diagrammatic NCA tech- 
nique [8, 11, 12]. We consider the case of infinite repul- 
sion, U — ¥ oo. It must be noted that different sets of 
self-consistent equations have to be solved in order to 
compute the equilibrium and non-equilibrium Green's 
functions. Therefore, computing the current from Eqs 
(2) and (3) allows us to test the validity of the methods 
used. 

For numerical evaluations, we consider a flat con- 
duction band with band width 2D and also take Vl = 
Vr = V, T u a = nV 2 /D. As it is clear from Eq. (2) 
and Eq. (3), the main dependence on the conductance 
at low enough temperatures T is given by the spectral 
weight close to the Fermi level. It is then useful to un- 
derstand the behavior of the spectral density of the QD 
for different conditions, specially in the Kondo regime, 
where an enhanced conductance is expected. In Fig.l 
we show the resulting p a for different applied bias volt- 
ages. We take Fl = 1 as our unit and set Ed = —4, 
D = 10 and T = 5 x 10" 5 . We fix the Fermi level 
ep — 0. Note that — Ed S> Fl so the localized level 
is always occupied ((n CT ) ~ 1). This corresponds to the 
Kondo regime, in which there is a localized spin inter- 
acting with the conduction electrons. 

At zero bias there are two peaks in the spectral den- 
sity. The one centered close to Ed, the charge transfer 
peak, is the result of a non-interacting orbital hybridiz- 
ing with a conduction band. If temperature were higher 
than the relevant low-energy scale of the problem, the 
Kondo temperature Tk, this would be the only peak 
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Fig. 3. Conductance as a function of temperature obtained 
from lineal response formalism (continuous curve) and differen- 
tial conductance curves evaluated at zero bias(dots) and NRG 
predictions. 

in the spectral density. Since T is very low, the low- 
energy physics is dominated by the Kondo singlet be- 
tween conduction electrons and the localized one. The 
localized spin leads electrons close to Fermi level to 
suffer spin-flip processes giving rise to a screening ef- 
fect. This is the reason for the increase of the spectral 
weight at €f shown in the figure, which corresponds to 
the Kondo peak. Its width is related with Tk- When 
bias is turned on (dashed curves in Fig. 1) the Fermi 
level of each metallic contact is shifted. We set /il = 
— fj.it = eV/2. This energy shift produces a splitting of 
the Kondo resonance since conduction electrons com- 
ing from both leads contribute to the screening pro- 
cess. As a direct consequence the spectral weight at the 
equilibrium Fermi level decreases and a lower conduc- 
tance is expected. In Fig. 2, we show the differential 
conductance for several temperatures. Since bias volt- 
age is applied symmetrically the current I(V) is an odd 
function and thus conductance dl /dV(V) an even one. 
dl/dV curves are Lorentzian-like with the maximum 
at zero bias. The value at the maximum and the width 
depend on temperature. 



Fig. 4. Conductance as a function of gate voltage V g for dif- 
ferent temperatures obtained from lineal response formalism 
(continuous curve) and from differential conductance curves 
evaluated at zero bias(rcd dots). 

From the maximum of the curves of Fig. 2, we can 
build a point-by-point curve of G vs. temperature. The 
result is shown in Fig. 3. We show also a continuous 
line curve which corresponds to an equilibrium calcu- 
lation of the conductance in the linear-response regime 
by means of Eq. 3. The dashed line curve is the em- 
pirical formula derived from the NRG calculations[13]. 
For high temperatures T 2> Tk, there is no Kondo res- 
onance and thus the spectral weight at the Fermi level 
is low. There is an intermediate region where thermal 
fluctuations compete, and at low enough temperatures 
T <C Tk, the Kondo effect is fully developed and con- 
ductance tends to a saturation value. However, at T <C 
Tk, the NCA overestimates the Friedel's sum rule and 
therefore the conductance exceeds the unitary limit. 

As it is shown in Fig. 3, there is an excellent agree- 
ment between the results from the non-equilibrium cal- 
culation and the equilibrium ones. Moreover, there is 
also a great correspondence to the results from NRG. 
We stress that the calculation of linear response con- 
ductance implies only equilibrium quantities while the 
out of equilibrium solving procedure is more complex 
and deals with lesser and greater Green functions. Since 
the first of this approaches is valid just under the con- 
dition of proportional couplings, the agreement we find 
is a useful check that supports reliability to the most 
general procedure based on the calculation of the cur- 
rent by means of Eq. (2) . 

We turn now the discussion to the conductance as a 
function of gate voltage V g . The energy of the localized 
level of the QD is proportional to this voltage eV g = 
—Ed and thus it is possible to perform a transistor- 
type experiment by the control of this parameter. In 
Fig. 4, we present the NCA results for different tem- 
peratures. The understanding of this outcome is di- 
rectly connected to our previous analysis of the spec- 
tral density. For high temperatures (black continuous 
line in the figure) conductance shows just a symmetric 
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peak centered at V g = 0. This corresponds to the lo- 
calized level placed at the Fermi energy, the optimum 
condition for conduction electrons to pass from left to 
right metallic contact. For temperatures lower than Tk 
(dashed curves in the figure) the behavior is completely 
different. As soon as the energy of the localized level is 
below the Fermi energy, the Kondo effect develops and 
the spectral density shows not only the charge transfer 
peak but also the Kondo resonance. This is the reason 
for the plateau of conductance as a function of gate 
voltage. Our results qualitative agree with those ob- 
tained previously with NRG [14] . We bring more details 
on the evolution with temperature. As it is shown in 
Fig. 4, at finite temperature conductance start to decay 
at some V g . This feature has to do with the fact that 
Tk ~ cxp{_Ed}, i.e. smaller for more negative values 
of Ed- Since temperature is finite, T turns bigger than 
Ti({Ed) at some point destroying the Kondo effect and 
the plateau. It must be noticed that for T < Tk, the 
NCA results are not reliable within the empty orbital 
regime, E d > 0, due to the appearance of a spike with 
non physical spectral weight at the Fermi energy. 

For T = 1.5 x 1CP 4 we show (with dots) in Fig. 4 the 
values of conductance obtained by the procedure stated 
previously from the calculations out of equilibrium. In 
the inset of the figure we show differential conductance 
curves for several V g . We observe that the maximum of 
the conductance keeps the same while the curves get 
narrower for greater values of gate voltage. This is a 
direct indication of the variation of Tk with V g . 

There are three different way to define the character- 
istic energy scale Tk from the physical magnitudes ad- 
dressed in this work: At equilibrium it can be obtained 
from the half-width at half of the maximum (HWHM) 
of spectral density for T — ► (T£). Out of equilibrium, 
from the HWHM of differential conductance curves, 
rpdi/dv p rom {.jjg equilibrium conductance, Tk is de- 
fined as the temperature such that G(T) = G(0)/2. 
From our results, T§ = 8.2 1CT 3 , T K = 9.1 1CT 3 = 
l.HTg and T K ' /dv = 12.1 1(T 3 = 1.48Tf. As ex- 
pected, the three values are of the same order of magni- 
tude. These realtions can be used to estimate any one 
of them in a case it might not be accessible and to test 
the model. 



4. Conclusions 

In this work we study the transport through an in- 
teracting quantum dot described with the Anderson 
model. We use the NCA to calculate the conductance 
as a function of bias and gate voltage, and temperature. 
We find a great agreement between the results from the 
non-equilibrium calculations and those from the linear 
response regime which imply only equilibrium quanti- 



ties. The results for conductance versus temperature 
also agrees with those from NRG calculations. We ana- 
lyze the conductance as a function of gate voltage and 
observe the formation of a plateau for low enough tem- 
peratures within the Kondo regime in agreement with 
previous results. As a consequence of finite tempera- 
ture, the conductance decays for a given value of the 
gate voltage destroying the plateau. We finally discuss 
several methods, in and out of equilibrium, which al- 
lows the determination of Kondo temperature. We find 
that the values obtained are of the same order and pro- 
vide numerical relations beetwen them. While the con- 
ductance at equilibrium of the model is known from 
NRG calculations, we show that the NCA is able to 
provide reliable results at a lower computational cost, 
and to extend the results out of equilibrium, for which 
few alternative techniques exist. Moreover, the NRG 
can miss structures in the spectral density which are 
not near to the Fermi energy [15,16]. This turned out 
to be important to explain a plateau observed in G(T) 
in Ceo QDs [2] on the triplet side of a quantum phase 
transition, using the NCA [16]. 
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